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Abstract. The standard phcnomcnological Hamiltonian of a small superconducting 
Josephson junction in the charge regime (Cooper Pair Box) produces a model of the 
effective charge qubit with possible applications to quantum information processing. 
In this note a new model based on the BCS Hamiltonian with individual tunneling 
yields an effective multi-level picture with a highly degenerated level placed between 
the ground state and the excited state. Unlike in the standard approach, the excited 
Cooper pairs play here an important role. For such a system coupled to a zero 
temperature bath the additional levels act as a probability sink. In contrast to 
the standard large-spin model the coupling to phonons can be an effective source 
of dissipation. This model provides also alternative explanations of various effects 
observed in experiments and sheds new light on the issue of Josephson junctions as 
macroscopic quantum systems. 
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In the last decade remarkable experiments were performed involving measurements 
and manipulations of states for a single or several nanoscopic Josephson junctions (JJ) 
which were consistently interpreted in terms of two-level quantum systems [1, 2]. The 
type of the J J considered here, the Cooper Pair Box (CPB), is a circuit consisting of a 
small superconducting island coupled via Josephson junction to a large superconducting 
reservoir [3, 4, 5, 6] . One of the standard forms of the CPB Hamiltonian is the Josephson 
Hamiltonian 

H = 4Ec{h-ng)^-§^J,. (1) 

Here, jk,k — x,y,z are j-spin operators, where j/2 is a an average total number 
of Cooper pairs on the island at equilibrium. The first term describes the Coulomb 
repulsion and the parameter Ug is an external control proportional to the gate voltage. 
The second term accounts for the tunneling of Cooper pairs between the island and 
the large electrode with the magnitude given by the Josephson energy Ej < Eq- At 
temperatures ksT << Ec, and putting < < 1, it is enough to restrict the analysis 
to the two lowest Coulomb energy states, say |m = 0), |m = 1), (for an integer spin) 
what leads to the picture of charge qubit. 

Although there exist derivation of the Hamiltonian (1) mainly based on the formal 
quantization of the phenomenological circuit equations or the two-mode approximation 
for the Bose-Einstein condensate(BEC) of Cooper pairs, the more rigorous approach is 
still needed. The main problem with the single degree of freedom model is its collective 
coupling to environment which should lead to a scmiclassical behavior with normal 
fluctuations ~ ^/J of the Cooper pairs number [8]. This effect is observed in atomic 
BEC systems [7] which arc described by the same Hamiltonian (1) but clearly not in 
CPBs where charge fluctuations are strongly suppressed. In this note a solution of this 
problem is proposed, which is based on the individual tunneling model combined with 
the individual coupling to a phonon's bath at zero temperature. The detailed derivations 
of the model Hamiltonians and the analysis of Markovian master equations, both for 
CPB and the current biased junction (phase qubit) will be presented in the forthcoming 
publication [12]. 

Reduced BCS Hamiltonian In the following a more fundamental approach is used, 
which is based of the simplifled BCS model [10, 11] which captures all essential features 
of superconductivity and should be also valid for small JJ's including CPB. 
In the absence of tunneling the island can be described by a system of Cooper pairs 
treated as hard-core bosons with the BCS mean-field interaction. Equivalently, one can 
consider a system of K spins-1/2 with spin operators s^, a — x,y,z satisfying 

[§1, sf] = iSkisl and cyclic permutations. (2) 

Defining the collective spin operators 

K 

J = {Jx,Jy,Jz), Ja = ^sl , a = x,y,z, (3) 

k=l 
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one can write the mean-field BCS Hamiltonian 

^-ci = -| E = -|(J^ - J'. + I). (4) 

k,l=l 

Here, = + is^, = sf — is^ create and annihilate a Cooper pair, K is an even 
number of electron levels in the cut-off region around the Fermi energy and g is an 
effective pairing potential. Using the decomposition of the Hilbert space for K spins- 
1/2 into subspaces corresponding to the irreducible representations of SU{2) of the 
dimension 2j + 1 and multiplicity 

K/2 

C2" = 0c2^+i<g,C'-^- (5) 

j=0 

one can use as eigenvectors of i^red the orthonormal basis \j, m; r) 

J^Ij, m; r) = j{j + m; r), J^|j, m; r) = m|j, m; r) (6) 

where m = + f — i,'^, One should notice that the collective 

spin operators (3) do not correspond to a single spin, like in the case of the Josephson 
Hamiltonian (1) but are direct sums of j-spin operators with j = K/2, K/2 — 1, ...,0, 
and with many copies for a given j < K/2. The states with j — K/2 — p describe 
excitations composed of p excited Cooper pairs. In contrast to unpaired electrons this 
type of excitations are not common in the literature. The main message of this paper 
is to argue that they play a crucial role in the description of CPB. 

As at zero temperature the energy levels up to Fermi level are filled the number 
of Copper pairs {K/2 + m) determined by the quantum number m is close to K/2 
and hence m ~ 0. For the finite temperature case one obtains a superconducting 
phase transition at the critical temperature ksTc — g/2 [11]. For a fixed m the energy 
difference between the unique ground state \K/2,m) and the {K — l)-degenerated first 
excited level \{K/2) — l,m;r) is equal to g what implies that g has a meaning of the 
zero temperature superconducting gap ((7 ~ A). 

Cooper Pair Box - BCS Hamiltonian 

Adding to the Hamiltonian ^red the term describing Coulomb repulsion one obtains 
the Hamiltonian of the isolated island 

H = - J,' + I) + 4Ec{l - mf (7) 

where the charging energy Ec fulfills ksT « Ec and fh ( tuq < m < tuq + 1, |mo| << 
K) determines the number of Cooper pairs in the systems, and can be controlled by an 
external voltage. For experimental realizations the parameters (7 ~ A ~ /c^Tc and Eq 
are of the same order of magnitude (~ 1 kelvin), while K ~ 10*^ — 10^. 

Due to the Coulomb blockade described by the second term in (7) one can consider 
only the states with m — tuq, tuq + 1 satisfying rh G [mo, mo + 1]. Taking into account 
that mo « K and extracting the irrelevant common constant 2Ec + 4£'c[(mo — rhY + 
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(mo — m)] — gK/A we have an effective Hilbert space spanned by tlie following vectors 
with the corresponding energies denoted by the simplified symbols 

|0) = \K/2, mo) Eo = -2Ec{l - 2ng), 



(8) 



|1) = \K/2,mo + 1) E,= 2Ec[l - 2ng - 4^], 

\s- 0) = \K/2 - 1, mo; s) = -2Ec[l - - 4~g], 

|r; 1) = \K/2 - 1, mo + 1; r) Wi = 2Ec{l - 2ng) 

where s,r — 1,2, K — 1 and 

rig^m-mo-g , g^ (9) 
Tunneling processes 

Tunneling process for Cooper pairs produces superpositions of states on the island 
which differ by a single pair and can be described by the reduced tunneling Hamiltonian 
acting on the states of the island 

1 ^ 

fre^-^J2((^kK + ^kK) (10) 
fc=l 

with, generally complex, tunneling probability amplitudes f3k- Introducing the total 
amplitude (3 — ^^=i one can decompose the tunneling Hamiltonian into collective 
and individual parts 

Tred = -(Re(/3)4 + Im(/3) J,) + -J^ [iPk - + {h - j^m] .(11) 

fc=i 

The collective part of (11) preserves the subspaces of a given j. Therefore, if the 
collective part dominates one could consider only the states with j — K/2 to obtain 
the standard picture of a large spin with the Hamiltonian (1). Indeed, the obtained 
collective Hamiltonian 

^ _ A(j2 _ j2 ^ ^ 4^^^ _ ^ i_(Re(/3)4 + Im(/3) 4) (12) 
K K 

restricted to the j = K/2,\m\ << K subspace and in the limit of large K is unitarily 

equivalent (up to the irrelevant constant) to the large spin Hamiltonian (1) with the 

Josephson energy equal to the collective component Ej — Ej = \f3\. To compare the 

magnitude of the collective component of the Josephson energy Ej with its individual 

counterpart given by 

- E 1/3. - - E - ^mr' (13) 

k=l k=l 

one can consider a simple toy model with (3k = Ae^'^'^, k — K — 1. Then 

M _ piXK\ 9141 , 

E'\^\A y\ < , ' L , E'^\A\^/K (14) 



what implies for a generic A that i?} ~ -E'jV K. 

On the other hand for purely random amplitudes (5^, = \Pk\'^ and therefore 
E'^j ~ Ey The real system should be placed between these two extremal cases of strong 
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interference and random behavior what imphes that the ratio Ej/Ej increases as a 
certain positive power of K leading to the domination of the individual coupling. 

For large JJ's with small Ec this effect is suppressed by the fact that the typical 
level splitting for a fixed j determined by the Coulomb repulsion is much smaller than 
the level splitting for different j's given by the superconducting gap. For small junctions 
those energy scales are comparable and the individual tunneling prevails. This implies 
that the matrix elements of Tred between the vectors with the same j are negligible in 
comparison with the elements between vectors with \j —f\ — 1. Hence the only relevant 
matrix elements are the following 

(r;l|ted|0) , (s;0|fred|l) =6- (15) 

Using (15) one can derive the full effective Hamiltonian of the CPB including (12,10) 
which is a direct sum of two similar terms 

HcpB = H^PB © Hl;pB (16) 

acting on the subspaccs 'H^g,'HlQ spanned by {|0),|r;l),r = 1,2, ...,K — 1} and 
{|1), |.s:0),s = 1,2,...,K — 1}, respectively. To a large extend both "subsystems" can 
be treated separately and completely analogically. Depending on the parameters cither 
T-L^Q or Hlff contains a ground state of Hcpb- The ground state is a starting point 
of the controlled evolution in all experiments and as shown below both subspaces are 
invariant with respect to the Hamiltonian and controls. Hence, in the following, only 
one subspace, say T-LeS^ represents the accessible system. 
Charge qubit picture Defining a normalized vector 

\0 = iEj)-'J^^r\r;l), Ej=(j^m'y" (17) 

r=l r=l 

one can introduce the vectors |±) 

l+)=cos^|0+sin^|0) , |-)=cos^|0)-sin^|0 , (18) 

where 9 is defined by cos^ = E{ng)/ ^/E(ji^yTEj, E{ng) = 4Ec(l - 2ng). The qubit 
observables are given by 

^+ = li^^'+^^n = , = , = i+)(+i+i-)(-i-(i9) 

Extracting from Hqpb overall constant E{ng)/2 + ^Ecn^g one obtains a new form 
of the Hamiltonian 

^Spb = + , ^ = sjE{ngf + El (20) 

with two eigenvectors |±) separated by the energy difference uj and the third (fT — 2)-fold 
degenerated level corresponding to the projector Pq == ^^^^ \r; l)(r; 1| — ^iid the 
energy E{ng)/2. This third level lies always between |ib), as — a;/2 < E{ng)/2 < uj/2. 

Notice that both parameters, Ej and g are related to the gap A and hence can 
be modified by the external magnetic field. The time-dependent external control is 
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performed by the coupling through the total electric charge operator Q and the total 
electric current J = i[HcpBj Q] which possess a similar direct sum structure as the 
Hamiltonian. When restricted to they read 

Q° = e[(sin - (cos e)a' - a° - 2Po] , J° = 2e a;(sin . (21) 

Obviously, if the system is completely isolated, then the qubit Hilbert space spanned by 
|±) is invariant with respect to the Hamiltonian and the external control yielding the 
usual model of charge qubit. The {K — 2)-fold degenerated energy level corresponding 
to Jo and the lowest energy levels of the second subsystem become important when the 
coupling to an environment is discussed. 
Dissipation and decoherence processes 

There exist many hypothesis concerning the leading mechanism of dissipation 
and decoherence in CPB. Some of the mechanisms can be, in principle, reduced by 
a proper engineering and the others seem to be not effective enough to account for 
the experimental data. The presented model of individual tunneling which couples the 
ground state to excited Cooper pairs states allows for a new mechanism of dissipation 
due to excited Cooper pair - phonon coupling. As the typical CPB frequency u lies 
well below the Dcbyc cut-off the dissipation through the coupling to phonons is not 
suppressed like, for example, in the case of quantum dots. The corresponding acoustic 
wavelength is of the order of 0.1 fim, comparable to the size of the CPB and to the 
average distance between electrons forming a Cooper pair, what implies rather individual 
coupling of different Cooper pair modes |A;) than the collective one. It follows again that 
the value of j is not conserved in the relaxation processes. The detailed analysis of the 
model based on the Markovian master equation for the reduced density matrix p of 
the CPB will be presented in [12]. Here, the written below Bloch equations for the 
level occupation probabilities p±. = (±|/)|±) , po = Tr(pfo) and the qubit coherence 
a — Tr{pa~^) can be treated as a phenomenological description: 

^^=-(r + rW + rWK, 

^ = T^iW - rfc'p., (22) 

^ = rp+ + T^wPo - 

^ =(i.-i[r + rw + rw + rt'i)a. 

The decay rates describe the relaxation and leakage processes depicted on the Fig.l. 
In the standard 2- level model only a single relaxation rate l/Ti corresponding to P in 
(22) is present, but the additional pure dephasing rate is added. Here pure dephasing 
is absent, because it cannot appear for a linear coupling to a bosonic bath (see the 
discussion in [14]) and the additional decay of the coherence amplitude a is due to the 
probability leakage. 

The Bloch equations (22) and their extensions including a resonant electromagnetic 
perturbation will be used in [12] to analyse the experimental data, both for charge and 
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Figure 1. Level structure and relaxation processes for CPB model, 
- qubit levels , highly degenerated level. 



phase qubits. Some observed effects like, for example, instability of the CPB ground 
state [4] and the two time scales of energy relaxation for phase qubits [13] will find an 
alternative explanations. 

Conclusions The idea that the collective tunneling of Cooper pairs is strongly 
suppressed by the interference, what for small Josephson junctions leads to the 
domination of an individual tunneling, yields a new model of a Cooper pair box. This 
model conceptually differs from the standard charge qubit picture but nevertheless 
reproduces experimental data and even provides simple explanations for some issues. It 
suggests also the existence of a new effective relaxation mechanism due to the interaction 
of excited Cooper pairs with phonons. The fundamental difference between the standard 
and the new model is that the later cannot be treated as a model of macroscopic 
quantum system with a single degree of freedom (e.g. large spin) but involves other 
degrees of freedom represented by highly degenerated states of excited Cooper pairs. 
This explains the absence of environmental effects leading to a semiclassical behavior 
for large quantum numbers. 
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